VICIOUS WALKERS AND RANDOM CONTRACTION 
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JONATHAN NOVAK 



Abstract. The ensemble CUE'''' of truncated random unitary matrices is a 
deformation of the usual Circular Unitary Ensemble depending on a discrete 
non- negative parameter q. cue'i' is an exactly solved model of random con- 
traction matrices originally introduced in the context of scattering theory. In 
this article, we exhibit a connection between CUE'''' and Fisher's random- 
turns vicious walker model from statistical mechanics. In particular, wc show 
that the moment generating function of the trace of a random matrix from 
cue'''' is a generating series for the partition function of Fisher's model, 
when the walkers are assumed to represent mutually attracting particles. 



1. Introduction 

1.1. Truncated random unitciry matrices. Fix an integer d > 1 , and let 
M(i(C) be the space of d x d complex matrices. Consider the linear map 



(1) 



T:: 



which acts by removing the last row and column of a matrix. For example, T 
M3 — > Ml is defined by 



(2) 



Zll 


Z12 


Z13 


Z21 


Z22 


Z23 


Z31 


Z32 


Z33 







Z12 




Z21 


Z22_ 



in operator norm for any 



T is called the truncation map. Since ||T(M)|| < 
M e Md+i , T maps the unitary group 

(3) Ud+i = {U e Md+i : U* = U-^ } 
into 

(4) B<i={PeMd:||P|| <!}, 

the semigroup of linear contractions of Euclidean space C^. More generally, for any 
integer q > 0, the q-fold composition T'''' maps each matrix in M^+q to its d x d 
principal submatrix and thus sends Ua+q into M^. 

The truncation operator induces a very natural deformation of the Circular Uni- 
tary Ensemble (CUE) from random matrix theory. Consider the unitary group 
U(j+q as a probability space, with the Borel cr-algebra and Haar probability mea- 
sure. The pushforward y^^^ of Haar measure on Ud-i-q under T'''' is a Borel prob- 
ability measure on B^, and thus one obtains a random matrix ensemble 

(5) {(Bd,Yli'"):d>l} 
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Figure 1. Three random-turns particles moving between ground states. 



which will be denoted CUE''^' and called an ensemble of truncated random unitary 
matrices. Since T'°' is the identity operator on M^, CUE'°' reduces to the usual 
Circular Unitary Ensemble of unitary matrices under Haar measure when q = 0. 

Ensembles of truncated random unitary matrices were first studied by Sommers 
and Zyczkowski |35j in the context of quantum chaotic scattering, and have been 
further investigated by Petz and Reffy [21], [31], Fyodorov and Khoruzhenko [Tf] . 
and Neretin |28| . They are an important technical ingredient in Krishnapur's recent 
study of random matrix- valued analytic functions p5'j . 

Certain averages over CUE are well-known to be related to the combinatorics 
of increasing subsequences in permutations |33j , [4] and vicious walkers on Z [T] , 
|15| . A natural question is whether (and in what form) this relationship with 
combinatorics extends to the deformed ensemble CUE'''' . 



1.2. Fisher's random-turns model. Vicious walker models were introduced in 
statistical mechanics by Fisher in order to model wetting and melting [T2]. Fisher's 
random-turns vicious walker model consists of a system of d particles ("walkers") 
initially occupying sites 

(6) ^ = m > M.2 > • • • > 

on the integer lattice Z. Note that the walkers are labelled from right to left, so 
that we list their positions in decreasing order. The system evolves in discrete time 
according to the following rule: at each instant, a single random particle makes a 
random unit jump left or right (a "random turn"), subject only to the constraint 
that no two particles can occupy the same lattice site simultaneously. The function 
Zd(N; p., A) which counts the number of ways in which d random-turns particles 
can depart initial sites \i and arrive at new sites 

(7) A = Ai > A2 > • • • > Ad 

at time N is the partition function of the model. 

If one assumes the existence of an attractive force between the particles, then 
the system can be at equilibrium only when the particles are on adjacent sites. In 
this case we use the simplified notation Zd(N; q) to denote the number of ways in 
which d random-turns particles can move between ground states q sites apart in 
N instants. Figure 1 gives an example of a sequence of configurations of mutually 
attracting random-turns particles counted by Z3(10;2). 
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1.3. Main result. Our main result is the following. Consider the matrix integral 

,xTr(P + P')-,,(q) 



H 



dx q 



(8) Gd(x;q) 
where 

i=0 

is the hook-product of the d x q rectangular Young diagram. Gd(x; q) is a scaled 
and shifted version of the moment generating function of Tr(p[j^' + Pd'''*), where 
pjj''' is a random matrix from CUE'''' . 

Theorem 1.1. Gd(x; q) is the exponential generating series for the partition func- 
tion of a system of d mutually attracting random-turns particles: 

(10) Gd(x;q) = X Zd(N;q)^. 

N>0 

Below, we use this interpretation of Gd(n-;q) as a generating series to give a 
completely combinatorial proof of the following identity, which expresses the matrix 
integral Gd(x; q) as a Toeplitz determinant of Bessel functions. 

Theorem 1.2. We have 

(11) Gd{x;q) = det(Iq+j_i(2x))i<i,j<d, 

where the entries in the determinant are modified Bessel functions. 

Remark: Neretin [28 has shown that, when q > d, the measure y[j''' becomes 
absolutely continuous with respect to Lebesgue measure on Bd and has density 

(12) .^^'^'^^ det(I-P*P)''-'^dP. 

Hdx(q-d) 

Therefore when q > d, Gd{x; q) can be written more concretely as 

gXTr(P + p-)det(I-P*P]q-ddP. 



(13) Gd(x;q) = ^j^ 



xdq 



dx(q-d) 

1.4. Connection with the increasing subsequence problem. Recall that a 
permutation cr from the symmetric group S(n) is said to have an increasing subse- 
quence of length k if there exist indices 

(14) 1 < 11 < ii < • • • < ik < n 
such that 

(15) o-{i,) < (T(i2) < ••• < cT{ik). 

Increasing subsequences in permutations were first studied by Erdos and Szekeres 
|11] in the 1930s in connection with a Ramsey-type problem for the permutation 
group. In the 1960s, Ulam IJUj raised the problem of determining the number 
Udln) of permutations in S(n) with increasing subsequence length bounded by 
d. This came to be known as the increasing subsequence problem. Stanley's ICM 
contribution |38j gives a comprehensive survey, from a combinatorial perspective, of 
the vast literature on increasing subsequences in permutations and related topics. 
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Forrester [TJ has observed the following connection between the increasing sub- 
sequence problem and Fisher's random-turns model: 



/2n\ 



(16) Zd(N;0) = 



)ud{Ti), if N — 2n for some n > 
0, otherwise 



We will present a new proof of (a more general version of) Forrester's result 
below. Our approach is based on the fact that the configuration space of Fisher's 
random-turns model naturally carries the structure of a graded graph with com- 
mutative raising and lowering operators. This approach was inspired by Stanley's 
notion of a "differential poset," defined as a graded graph which arises as the Hasse 
graph of a poset and whose raising and lowering operators satisfy the Heisenberg 
commutation relation |36j . In particular, this observation implies that specializing 
q = in Theorems |1.1| and |1.2| we obtain 

(17) Gd(x;0) = y Ud(n)4^ = det(Ij_i{2x))i<t,j<d. 

^ — n!n! - 

n>0 



The identity between the matrix integral and the generating series of Ud(n.) is due 
to Rains [33], while the identity between the series and the determinant is due to 
Gessel [TH] . Both were discovered in the context of the enumeration of permutations 
with bounded increasing subsequence length. 



2. Configuration Space as a Graded Graph 

2.1. Configuration space of random-turns particles. The configuration space 
of a physical system is the set of all possible positions of its constituents. In the 
case of the random-turns model, this is the discrete space 

(18) >Vd={{Ai,A2,...,Ad)eZ'^:Ai > A2 > • • • > Ad}. 

Wd is familiar as a type A Weyl lattice, i.e. it is the intersection of Z'^ with an 
open Weyl chamber for the type A root system in M'^ (see e.g. [8] about root 
systems and their Weyl chambers). Wd becomes a simple, connected, locally finite 
graph when we declare vertices A e Wd adjacent if and only if they are unit 
Euclidean distance apart. The partition function Zd(N; p., A) of the random-turns 
model counts the number of walks of length IM from |j, to A on the graph Wd- 

Wd is also an example of a graded graph. For our purposes, it suffices to define a 
graded graph to be a pair (G, r) consisting of a simple, connected, locally finite graph 
G together with a function r : G — > Z, called the rank function, which associates an 
integer to each vertex of G. The incidence relation and the rank function on G are 
required to be compatible in the sense that if u,v G G are adjacent vertices, then 
either r(v) = r(u) -|- 1 (denoted u ^ v) or r(v) = r(u) — 1 (denoted u \ v). Graded 
graphs occur frequently in combinatorics as the Hasse graphs of posets (see [36] 
and '13]) and in representation theory where they are known as branching graphs 
or Bratelli diagrams (see e.g. [24 ). 

If we define a rank function r on Wd by 

d 

(19) r(A) =r(Ai,...,Ad) 

i=l 




RANDOM CONTRACTIONS 



5 



then (VV(i,r) is a graded graph. Note that we have normalized the rank function 
in such a way that the rank of the canonical ground state p = (d, d — 1 , . . . , 1 ) in 
Fisher's model is 0. 

Let (G,r) be a graded graph. By the unrefined partition function of G we mean 
the number Zg(N; u, v) of walks on G from u to v. We will refer to the number 
Zg (L^^R'^" . . . L^i R'^' ;u, v) of walks on G from u to v of the form 

(20) u // - --/\\ - --\ . . . V 

ai bi aic bic 

as the refined partition function of G. In the special case G = Wd we will continue 
to use the notation Zd(N;(a., A) = Zwj(N;|a, A) from the Introduction, as well as 
Zd(L'"<R'"' ...L^'R-^'-.^.A) := ZwJL^-'R'"' . . . L^' R"^' ; ^i, A). 

L and R as defined above are simply inert symbols recording instances of positive 
and negative steps on the configuration graph (cquivalcntly left and right jumps of 
the particles). We may represent L and R as linear operators acting on C[G], the 
free C- vector space spanned by the vertices of G, by defining 

(21) R(u) ^ V and L(u) := ^ w 

and extending linearly over C[G]. R and L are then called the raising and lowering 
operators on G. The monoid {L, R}* generated by the raising and lowering operators 
determines the refined partition function, since 

(22) ZgIL^^R"" ...L^^R"^u,v) = ML^^R"" . . . L^^ R"Uu) 
where [v] is the "coefiicient of v" functional on C[G]. 

Theorem 2.1. The raising and lowering operators on (Wd,T) commute. 

Proof. Let Wd be the graded graph with vertex set 

(23) Wd={(A,,A2,...,Ad) eZ^iAi > A2 > • • • > Ad), 
incidence relation as with Wd, and rank function defined by 

r 

(24) f(A)=f(Ai,A2,...,Ad):=^Ai. 

i=1 

Since the translation 

(25) A A := A - p 

is a rank-preserving isomorphism of Wd with Wd, it sufiices to prove that the 

raising and lowering operators associated to {Wd,f) commute. 

Let ^, A G Wd be arbitrary vertices. Consider separately the cases H- ^ A and 
H = A. 

Ccise (J. ^ A: [A]LR(|a.) counts the number of walks on Wd from (x to A of the form 

(26) 10. [0. + \ n- + Ci - Cj = A, 

where ei and are standard basis vectors of M.^. If such an i and j exist, 
then i ^ j by the assumption |J. 7^ A. Thus 

(27) H \ |i. - ej ^ - ej + et 

is a valid walk on Wd from \i to A, and we have a bijection 

(28) ^^y\-\ ^ ^\y' 
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between up-down walks on Wa from |j. to A and down-up walks on Wa fr'om 
]x to A. Hence [A]LR(m,) = [A]RL(h). 
Case |x = A: In this case, let r be the number of distinct entries of |j. = A, i.e. 
ii < i2 < • • ■ < It satisfy |J.i, > [Xi^ > ■ ■ ■ > [li^. Then it is clear that 
[^]RL(^] = [h]LR(h) =r. 

□ 



2.2. Determinant identities. Using Theorem |2.1| we can easily prove determi- 
nantal identities for the refined partition function of random-turns particles. Let 
Ik be the modified Bessel function of order k, 

(29) Ik(2x) -.^ZruTZ 



n>0 



r(n+l) r(n + k+1)' 



It is known that I_k(2x) — Ik(2x) for k G Z (see e.g. about Bessel functions 
and their properties). 

Proposition 2.2. Let ^, A S Wa be configurations with r[\i] < t(A), and let 
^1 > ■ ■ ■ i^n, ■ • • €E {L, R}* be any sequence of words in the raising and lowering 
operators on Wa verifying 

(30) degL = n and dcg^ = n-|-r(A)-r(|j.). 
Then 

^2n+T(A)-T(H) 

(31) ^Zd(W^;^,A) ^ =det(lA,_^,(2x))i<t,j<d. 

Proof. Since Wd is the intersection of Z'^ with an open type A Weyl chamber in 
K*^, it follows immediately from the Andre-Gessel-Zcilberger reflection principle 
[20] . [21], [19] that the generating series for the unrefined partition function is 

(32) Y. Zd(N;H,A)— =det(lA,_^J2x))i<t,j<d. 

N>0 

Now, an N-step walk from |j, to A on Wa exists if and only if N = 2n -|- r(A) — r(|j,) 
for some n > (this n being the number of negative steps). By Theorem 2.1 the 
number of such walks is 

(33) I ^ IZd(Wn;H,A) 

for any € {1, Rj* with degL = Ti and deg^ Wn + r(A) —r([i). Thus 
(34) 

> Zd(N;M.,A)— = > Zd(W^;H,A)— — — — — — 

^ — N ^ — V n / 2n-Fr A — r u. 

N>0 Ti>0 ^ ^ ^ I I J ^.I^IJ 

2T^+T(A)-r(^) 

(35) =Z^d(W,;kL,A)— — — — 

n!(n r{A) -r(^))! 

and the result follows. □ 



Corollary 2.2 generalizes recent results of Xin (i42j. Theorem 11 and Proposition 
12), who proves the case = L"^R^ using the "Stanton-Stembridge trick." 
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2.3. Young tableaux and increasing subsequences. Let Y be the Young graph, 
i.e. the Hasse graph of the lattice of Young diagrams partially ordered by inclusion 
of diagrams (see f57J). Y is a graded graph, where the rank |A| of a Young diagram 
A e Y is the number of cells in A. Y has a unique vertex of rank (the "empty 
diagram"), and all other vertices have positive rank. Walks on the Young graph 
are known as oscillating Young tableaux. 

Let be the induced subgraph of Y whose vertices are the Young diagrams 
with at most d rows. It is a well-known consequence of the RSK correspondence 
(see [37] as well as the original article of Schensted |31J) that 

(36) ZyJL'^R^-JJ) =Ud{n), 

where Uafn.) is the number of permutations in S(n) with no increasing subsequence 
of length greater than d. 

Observe that there is a canonical embedding i : Y^ w obtained by mapping 
each Young diagram A in Y^ onto the vector of its row lengths. Composing i with 
the translation 

(37) A A° A + p 

gives an embedding of Yd into Wa- Thus it is an immediate consequence of Theorem 
Othat 

^(2-)ud(n), if N =2n 
0, otherwise 



(38) Zd(N;0) = 



which is precisely Forrester's result (16 1. We therefore have the following Proposi- 



tion as a direct consequence of Proposition 2.2 
Proposition 2.3. For any d > 1 , 



^2n 



(39) V Ud(n)— =det(Ii_j(2x))i<t,j< 

— n In I ' — — 



T1>0 



Proof. Choose |x = A = p in Corollary |2.2[ □ 



This is precisely Gessel's identity, the original proof of which appears in [18] . 
Alternative proof of this result were later given by Gessel, Weinstein, and Wilf 
^41j , Tracy and Widom [33], and Xin [35]. Gessel's identity was the starting point 
of Baik, Deift, and Johansson ^ in their groundbreaking work on the limiting 
distribution of the length of the longest increasing subsequence in a large random 
permutation. 



3. Proof of the main theorem 



We now give the proof of Theorem |1.1[ which links truncated random unitary 
matrices with Fisher's random-turns model. The proof is based on the integral 
identity 



(40) 



Gd(x;q) = 



,xTr(U + U* 



det(U*)''dU 



which we deduce from the following remarkable result of Wei and Wettig. 
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Theorem 3.1. ^ For any matrices X, Y e Mia+qjxd verifying det(Y*X) ^ the 
following holds: 



(41) 



,Tr(Y'UX+X*U*Y) 



dU = Hd> 



To get (40 1 from Theorem 3.1 choose 

ft 

t 



(42) X=Y= 



-Ti(UY*Y+X'XU*) 



det(UY*X)-'^dU. 



_0 ... 

where t is a non-zero real number. Then Y*X — t^I, and Y*UX = T''''(U) for any 



unitary matrix U G U^+q. Thus Theorem 3.1 reduces to 

Tr(P + P*).,,(q) 



(43) 



t2dq 



H 



dx q 



'y:i"'(dP) = Gd(t^;q) = 



Tr(U+U' 



det(U*)''dU, 



vahd for all real t 7^ 0. This identity also holds true at t = 0, since then the left 
hand side is obviously equal to 0, and the right hand side is also equal for the 
following reason: 



(44) 



F(U)G{U*)dU = 



for any two homogeneous polynomials F, G in the entries of U and U* with deg F ^ 
deg G (this is a standard property of Haar measure, see e.g. [9] for a proof). Thus 
(43 1 is true for all real numbers t and we obtain (40). 



To complete the proof of Theorem |1.1| we compute the power series expansion of 
the unitary matrix integral ( 40 ) using techniques from symmetric function theory 



(see |37j, Chapter 7). Recall that a monotone walk 
(45) A 



on the Young graph Y is called a standard Young tableau of shape A. Following [37], 
we will use the notation 



(46) 



f 



A 



Zy(RI^I;0,A) 



for the number of standard Young tableaux of shape A. 

Expanding the integral (|40|) as a power series in x, we have 



(47) 



Gdl'Ciq) 



m,n>0 



(TrU)"^(TrU*)^det(U*)''dU- 



mln! 



Appealing again to (|44|), this becomes 
(48) 



Gd(x;q) = Id{rL;q) 



^2n.+ dq 



n>0 



where 
(49) 



Id(n.;q) 



n!(n+ dq)!' 
(TrU)^+'^''(TrU*)^det(U*]''dU. 
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Now we evaluate the integral Iain; q). Let Ad be the C-algebra of symmetric 
polynomials in d indeterminates, and define the Hall scalar product on by 



(50) (f|g) 



f(U)g(U*)dU 



where f (U), g(U* ) denote symmetric polynomials f , g G evaluated on the spectra 
of U and U* respectively. The Hall product is a symmetric, non-negative bilinear 
form on A^ (see [37]). It is well-known that the Schur polynomials {sAlAGYd consti- 
tute a linear basis of A^, and moreover since each map U i— > Sa(U) is the character 
of an irreducible polynomial representation of Ud we have 

(51) (sa|Sh)=6a,h 

by Schur orthogonality (see e.g. [8]). 
Now 

(52) Id(n;q) = (e;^+'^''|e7e^), 

where ei , Cd G Ad are the first and last elementary symmetric polynomials. In 
terms of Schur functions, on has the linear expansion 



(53) = Z ^^'^ 



|A|=T1 

for powers of the first elementary symmetric polynomial (see [37]). It is also clear 
from the combinatorial definition of Schur polynomials in terms of semistandard 
Young tableaux that — Sdxq- Hence we have 



d — •'dx q ■ 

(54) Id(n;q)= ^ f^f^s^|sASdxq/ 

H.|=n+dq,|A|=n. 

The following property of the Hall scalar product is well-known: 

(55) (s^|sASv) = (s^/t,|sa), 

where s^^/^ is a skew Schur polynomial, which is by definition the zero polynomial 
unless |j, D "V. In other words, the adjoint of the "multiplication by s-v" operator is 
the "deletion of V operator. Hence we have 

(56) Id(n;q)= ^ f^f'(s^/dxq|sA) = ^ ^'^^''^^^ 

H.AeYd AeYd 
Vi|=n+dq,|A|=n. lAl— n 

where A -|- d x q is the concatenation of the Young diagram A with the rectangular 
diagram d x q. Thus 



(57) Id(n;q]= ^ Zy, (R-+'i'';0, A + d x qjZy, (R-;0, A) 

AGYd 
|A|=n 

(58) = ^ ZY,(R^+'^^;0,A + dx q)ZY,(L^;A + dx q,dx q) 



AGYd 

|A|=Tl 



(59) =ZYd(L'^R"+'''';0,dx q) 

(60) 'Zd(2n+dq;q), 
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and we conclude that 



(61) 
(62) 
(63) 
(64) 



n>C 



,2n+dq 



n>0 



In + dq 
n 



n+ dq)!n! 
1 



Y_ Zd(2n+ dq;q; 



n>0 



Zd(2n+ dq;q) 

-^2n+dq 

{2n+ dq)! 



^2n + dq 



(n+ dq)!n! 



= Y_ Zd(N;q 



N>0 



N! 

4. ASYMPTOTICS 



Let us now extract an explicit formula for Z(i(N; q) in terms of random contrac- 
tions from Theorem 11.11 We have 

^dq 

(65) 



Gd(x;q) = 



H 



dx q 



gXTr(P + p-)Y(^P)(dP) 



(66) 
(67) 

(68) 

Thus 
(69) 

Zd{N;q) 



N>0 

z 

n>0 

z 

n>0 



,dq 



H 



dx q 



(TrlP + P^jj^Y^i^'ldP) — 



H 



dx q 



|XrP|2n^(q) 



^2n+dq 



:dPi 



(2n)! 



-, by (441 



( 2n+ dq) 

Idxql 



ITrPI^-yfj^'ldP) 



^2n+dq 



(2n+dq) 







N! 





C^lBjTrP|^>:i^'(dP), ifN=2n- 
0, otherwise 



dq 



From this expression, we can easily obtain an asymptotic form for Zd(2n + dq; q) 
in the large q limit with d, n fixed. This corresponds to the "walk to infinity" with 
bounded backtracking Fisher's model. 

Recall the following classical result of E. Borel [5]. Let p — 
(pi , . . . ,Pd,Pd+i > ■ • • )Pd+q) be a uniformly random point from the real sphere 
gd+q-i ^ K'i+q_ xhen, as q — ) oo, ^qpi , . . . , y^Pd converge weakly to an inde- 
pendent family of standard Gaussian random variables. There is an ample gener- 
alization of this result to the setting of truncated random unitary matrices, first 
established by Petz and RefFy in [31] . 

Theorem 4.1. Let Pjj''' be a random matrix from CUE'''' . q — > oo, the reseated 

random matrix y'qPd''' converges to the random matrix whose entries are i.i.d. 
standard complex Gaussians (in the sense of pointwise convergence of spectral cor- 
relation functions). 

We remark that this result of Petz and Reffy follows readily from the work of 
Sommers and Zyczkowski [3S]. In particular, one knows from [3S] the the spectrum 
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of a random matrix P^j^' is a determinantal point process in the unit disc D = {z e 
C : |z| < 1 } governed by the Sommers-Zyczkowsi kernel 

(70) SZ<^'(z,w) = ^(1 - Izh^d - Iwh^ X ('^ ^0"'-'. 

^ j=o ^ ^ ' 

for any q > 1 . One then observes that the re-scaled kernel q^^ SZlj''' (q^^z, q^iw) 
converges to the Ginibre kernel 

I _ 

-z'w' 



1 "^"^ 1 

(71) Gina(z,w) = -e-<l^l'+l-l'' V - 

7T ' 1 I 



7T l! 

j=0 ' 



as q — > 00. An elementary geometric proof of Theorem 4.1 can be found in [26] . 
while a stronger assertion has recently been proved by Krishnapur j25j . 
Theorem |4. 1| immediately implies that 

(72) [ |TrP|^>l,^'(dP)^ ''^l"^+-;+"'^l'"^ 

as q — > oo, where zi , . . . ,Zrt are independent standard complex Gaussians and E 
denotes expected value. Since the moment sequence of a standard complex Gaussian 
z is well-known to be 

(73) E[z™zT =6^,nn!, 

and zi -|- • • • + za is a complex Gaussian of variance d, we have 

(74) 



|TrP|2-Y'^'(dP) 



as q — > oo. Applying this asymptotic form in (69), it follows from Stirling's approx- 
imation that 

(75) Zd(N; q) ^ (Ztt)^ ( ]]( il) d^n+d. + i ^n+i^ 

as q — > oo, where N — 2n + dq and n > 0, d > 1 are fixed but arbitrary. This 
asymptotic form reveals an interesting phase transition in the large q asymptotic 
behaviour of Z^fN; q) at n « 2^ from 

exponential in q 

76 — ^ 

polynomial m q 

to 

(77) (exponential in q) (polynomial in q). 
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